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Abstract 

Within the Hadron-String Dynamics (HSD) transport approach we study the sup- 
pression pattern of charmonia at RHIC with respect to centrality and rapidity em- 
ploying various model concepts such as variants of the 'comover absorption' model or 
the 'charmonium melting' scenario. We find that especially the ratio of the forward 
to mid-rapidity nuclear modification factors of J/^ {R^^l'"'''"^{J /^) / Rf^{J /^)) 
cannot be explained by the interactions with 'formed' comoving mesons or by the 
'color screening mechanism' alone. Only when incorporating interactions of the c or 
c quark with a pre-hadronic medium satisfactory results are obtained. A detailed 
comparison to the PHENIX data demonstrates that non-hadronic interactions are 
mandatory to describe the narrowing of the J/^' rapidity distribution from pp to 
central Au-|-Au collisions. The ^' to J/^' ratio is found to be crucial in disentan- 
gling the different charmonium absorption scenarios especially in the RHIC energy 
range. Furthermore, a comparison of the transport calculations to the statistical 
model of Gorenstein and Gazdzicki as well as the statistical hadronization model 
of Andronic et al. shows differences in the energy dependence as well as centrality 
dependence of the J/'^ to pion ratio which may be exploited experimentally to dis- 
entangle different concepts. We find additionally that the collective flow of charm 
in the HSD transport appears compatible with the data at SPS energies but sub- 
stantially underestimates the data at top RHIC energies such that the large elliptic 
flow V2 of charm seen experimentally has to be attributed to early interactions of 
non-hadronic degrees of freedom. 
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1 Introduction 



An investigation of the formation and suppression dynamics of J/\l/, Xc and 
mesons opens the possibility to address fundamental questions about the 
properties of the state of matter at high temperature and density. Up to date, 
a simultaneous description of the seemingly energy-independent suppression of 
J/\E' together with its narrow rapidity distribution and a strong elliptic flow V2 
of charmed hadrons - as found at the Relativistic-Heavy-Ion-CoUider (RHIC) 
- has presented a challenge to microscopic theories. The large discrepancies 
of present studies are striking in view of the success of the hadron-string 
transport theories in describing charmonium data at SPS energies. This has 
lead to the conjecture that the sizeable difference between the measured yields 
and transport predictions is due to a neglect of the transition from hadronic to 
partonic matter, e.g. a strongly-coupled Quark-Gluon-Plasma (sQGP). In the 
present work, we report new results on the charmonium nuclear modification 
factor RaAi rapidity distribution, the elliptic flow V2 of D mesons, the ratios 
{J/'^)/{tt) and '^'/{J/'^) for energies from about 20 A-GeV - relevant for 
the future Facility-for-Antiproton-and- Ion-Research (FAIR) - up to top RHIC 
energies. 

We recall that in the early stage of the nucleus-nucleus collisions the disso- 
ciation and the regeneration of J/\l/ by fundamentally different mechanisms 
are possible: The cc pairs produced early in the reaction - by gluon-gluon fu- 
sion in primary nucleon-nucleon interactions - might be completely dissociated 
in the dense medium and not be formed as bound states due to color screening. 
In this model scenario charmonia have to be recreated by some mechanism to 
yield a finite production cross section of J /'^ and ^' . The cc pairs might also 
be formed in some pre-hadronic resonance (color-dipole) state that will further 
develop to the charmonium eigenstates in vacuum. Such resonance states can 
be dissociated in the medium due to interactions with other degrees of freedom 
but also be recreated by the inverse reaction channels. Independently, char- 
monia might also be generated in a statistical fashion at the phase boundary 
between the QGP and an interacting hadron gas such that their abundance 
would appear in statistical (chemical) equilibrium with the light and strange 
hadrons PPI2]- In the latter model the charmonium spectra carry no informa- 
tion on a possible preceeding partonic phase. Indeed, in Ref. |3] a success of 
the statistical hadronization model [HSJ has been put forward. Another alter- 
native is the model for coalescence of charmonium in the sQGP For further 
variants or model concepts for charmonium suppression/enhancement we refer 
the reader to the reviews [7115] . In this work our aim is to shed some light on 
various model concepts by exploiting relativistic microscopic transport theory. 

The Hadron- String-Dynamics (HSD) approach [9J provides the space-time ge- 
ometry of nucleus-nucleus reactions and a rather reliable estimate for the lo- 
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cal energy densities achieved, since the production of secondary particles with 
hght and single strange quarks/ ant quarks is described well from SIS to RHIC 
energies [lOJ. As we will show in Section 2, the high energy-densities reached 
in Au + Au collisions at RHIC clearly indicate that a strongly interacting QGP 
(sQGP) has been created for a couple of fm/c in the central overlap volume. 
However, a proper understanding of the transport properties of the partonic 
phase is still lacking. Presently the effective degrees of freedom in the sQGP 
are much debated and charmonia are a unique and promising probe that is 
sensitive to the properties of the early (and so far unknown) medium. 

In the present systematic study we first test the HSD results for charmonium 
production in p + p and d + Au reactions at RHIC energies in comparison 
to the recent data. This is crucial in order to obtain an accurate baseline for 
the study of any anomalous suppression of charmonia in nucleus-nucleus col- 
lisions (see Sections E] and H]). The interactions of J/\I''s with mesons in the 
late stages of the collision (when the energy density falls below a critical value 
of about 1 GeV/fm^ corresponding roughly to the critical energy density for a 
parton/hadron phase transition) gives a sizable contribution to its anomalous 
suppression at all beam energies as demonstrated in Refs. [TT|12lll3lll4|15j . 
Accordingly, this more obvious 'hadronic' contribution has to be incorpo- 
rated when comparing possible models for QGP-induced charmonium sup- 
pression to experimental data. On the other hand, as known from our studies 
in Refs. [TTp2] charmonium interactions with the purely hadronic medium 
alone (which is modeled rather precisely by HSD) are not sufficient to de- 
scribe the J/\Ef suppression pattern at RHIC in detail. 

Based on the microscopic HSD transport theory, we investigate in particular 
the following scenarios for the anomalous absorption of charmonia: 

(1) the 'threshold melting' mechanism; 

(2) a dissociation by the scattering on hadron-like correlators, i.e. the 'co- 
mover' scenario; 

(3) additional scattering of charm with pre-hadrons which might be consid- 
ered as color neutral precursors of hadronic states (cf. Refs. p!6|17fl8|[T9] ). 




All implemented scenarios will be described in detail in Section O In Section [5] 
we will investigate in particular the effect of the interactions of charm quarks 
in the pre-hadronic medium on RAA{y) of J/"^ by comparing our calculations 
to RHIC data. To complete our study, we will provide excitation functions for 
the J/^ survival probability ^(J/^) and the ratios B^^a{J /^!) / B'^^a{^!') in 
Section [71 Furthermore, by studying the J/'^ to ti ratio as a function of the 
number of participating nucleons Npart, we will test the assumption of char- 
monium production by statistical hadronization as advocated in Refs. [2][3ll20] 




(subsection l7.2p . A summary of results as well as a discussion of open problems 
will close our study in Section 8. 
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Fig. 1. The energy density e{x = 0, y = 0,z;t) from HSD for a central Au+Au 
collision at ^/s = 200 GeV. The time t is given in the nucleon-nucleon center-of-mass 
system. 

2 Energy density 



The HSD transport model - employed in a large variety of ir+A, p+A, d+A and 
A+A reactions - allows to calculate the energy- momentum tensor T^y{x) for all 
space-time points x and thus the energy density e{x) = Too{x) in the local rest 
frame. In order to exclude contributions to T^jy from noninteracting nucleons 
in the initial phase all nucleons without prior interactions are discarded in the 
rapidity intervals [7/^^-0.4, 7/^^-^0.4] and [ypro-0.4:,ypro + 0A] where 7/^^ and 
i/pro denote projectile and target rapidity, respectively. Note that the initial 
rapidity distributions of projectile and target nucleons are smeared out due 
to Fermi motion by about ±0.4. Some comments on the choice of the grid in 
space-time are in order here: In the actual calculation (for Au+Au collisions) 
the initial grid has a dimension of 1 fm x 1 fm x 1/ 'jcm fm, where 7cm denotes 
the Lorentz 7-factor in the nucleon-nucleon center-of-mass system. After the 
time of maximum overlap tm of the nuclei the grid-size in beam direction 
Azq = l/jcm [fm] is increased linearly in time as Az = Azq + a{t — tm), where 
the parameter a is chosen in a way to keep the particle number in the local 
cells of volume AV{t) = AxAyAz{t) roughly constant during the longitudinal 
expansion of the system. In this way local fluctuations of the energy density 
e{x) due to fluctuations in the particle number are kept low. Furthermore, 
the time-step is taken as At = 0.2Az{t) and increases in time in analogy to 
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Az{t). This choice provides a high resolution in space and time for the initial 
phase and keeps track of the relevant dynamics throughout the entire collision 
history. 

The energy density e{r; t) - which is identified with the matrix element Too(r; t) 
of the energy momentum tensor in the local rest frame at space-time (r, t) - 
becomes very high in a central Au+Au collision at ^/s = 200 GeV as shown 
in Fig{T] (c/. Fig. 1 of [11] for the corresponding energy density evolution in 
case of central collisions at top SPS energies). Fig. [1] shows the space-time 
evolution of the energy density e{x = 0,y = 0, z; t) for a Au+Au collision at 
21300 AGeV or ^/s = 200 GeV. It is clearly seen that energy densities above 
16 GeV/fm^ are reached in the early overlap phase of the reaction and that 
e{x) drops after about 6 fm/c (starting from contact) below 1 GeV/fm^ in 
the center of the grid. On the other hand the energy density in the region of 
the leading particles - moving almost with the velocity of light - stays above 
1 GeV/fm^ due to Lorentz time dilatation since the time t in the transport 
calculation is measured in the nucleon-nucleon center-of-mass system. Note 
that in the local rest frame of the leading particles the eigentime r is roughly 
given by r ^ t/lcm- As seen from Fig. 1, the energy density in the local rest 
frame is a rapidly changing function of time in nucleus-nucleus collisions. For 
orientation let us recall the relevant time scales (in the cms reference frame): 

- The cc formation time ~ 1/M_l is about 0.05 fm/c for a transverse mass 
of 4 GeV; the transient time for a central Au+Au collision at y/s = 200 GeV 
is tr ~ 2RA/'~icm ~ 0.13 fm/c. According to standard assumptions, the cc pairs 
are produced in the initial hard NN collisions dominantly by gluon fusion in 
the time period t^. In fact, the formation time is significantly smaller than 

which implies the c or c quarks may interact with the impinging nucleons 
of the projectile or target for times t < t^. 

- Using the Borken estimate for the energy density and employing the time- 
scale tr = 0.13 fm/c, the energy density - after the nuclei have punched 
through each other - amounts to about 5/0.13 > 30 GeV/fm^ as quoted also 
in the HSD calculations in Ref. [H]. Even when adding the cc formation time, 
this gives an energy density ~ 5/0.18 ~ 28 GeV/fm^. So these numbers agree 
with transparent and simple estimates (cf. Fig. [1]) and illustrate the high initial 
densities after cc production from primary interactions. 

The energy densities quoted above are considerably different from the Bjorken 
estimate 

^7-1 ^ dN 



where < Et > is the average transverse energy per particle, dN/ drj the number 
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Fig. 2. Left part: The energy density per pseudorapidity interval dr] divided 
by the number of participant pairs (O.bNpart) from HSD (sohd hne) in comparison 
to the PHENIX data (dots) [21]. Right part: The Bjorken energy density ebj ■ t 
from HSD (sohd hne) for Au+Au cohisions at ^/s = 200 GeV in comparison to the 
PHENIX data (dots) [H]. 

of particles per unit of pseudorapidity, and r a formation time parameter 
often used as r = 1 fm/c. Furthermore, ttE^ denotes the overlap area for the 
corresponding centrality. Is is important to point out that the estimate ([1]) is 
only well defined for the product Te^j! The question naturally arises, if the 
transport calculations follow the corresponding experimental constraints. 

To this aim we show dET/drj (divided by half the number of participants Npart) 
from HSdH (l.h.s.) in comparison to the measurements by PHENIX [21]. Ac- 
cordingly, the Bjorken energy density esj - multiplied by the time-scale r ([1])- 
from HSD is shown additionally in the r.h.s. in comparison to the PHENIX 
measurements as a function of Npart- The similarity between the calculated 
quantities and the experimental data demonstrates that the space-time evo- 
lution of the energy-momentum tensor T^^^, in HSD is sufficiently well under 
control. We now may step on with the actual investigation of the charmonium 
dynamics. 



3 Implementation of charmonium production 



In order to examine the dynamics of open charm and charmonium degrees 
of freedom during the formation and expansion phase of the highly excited 
system created in a relativistic nucleus-nucleus collision, one has to know 
the number of initially produced particles with c or c quarks, i.e. D,D,D*, 
D\ Ds, 5„ D*, D*, J/^llS), <i/'{2S), XcilP)- In this work we follow the previ- 

^ The open source code is available from Ref. [22] 



6 




Fig. 3. Cross section for the differential J/^' production in rapidity (times the 
branching ratio to di-muons -B/^^) in pp colhsions at ^/s = 200 GeV. The HSD 
(input) parametrization (sohd hne) is compared to the PHENIX data (symbols) 
from Ref. |23j . 

ous studies in Refs. P)llfl3fl4|2^ and fit the total charmonium cross sections 
= Xc, J/^ 1 from NN collisions as a function of the invariant energy y/s 
by the expression [12] 



or{s) = ha il-^X f^V^(v^-v^), (2) 



S J yiT^i , 

where rrii denotes the mass of charmonium i while ^/soi = rrii + 2mN is the 
threshold in vacuum. The parameters in ([2]) have been fixed to describe the 
J/^ and data up to the RHIC energy = 200 GeV (c/. [II]). We use 
a = 0.16 mb, a = 10, j3 = 0.775. Note, that for the present study we updated 
our parametrization (i.e. reduced the parameter a by ~ 20%) according to the 
latest PHENIX data 



The parameters fi are fixed as f^^ = 0.636, fj/q, = 0.581, /^,' = 0.21 in order 
to reproduce the experimental ratio 

measured in pp and vrA^ reactions [25)126] as well as the averaged pp and pA 
ratio (-B^^(\l/')cT*')/(i?^^(J/\l/)(jj/^) ~ 0.0165 (c/. the compilation of exper- 
imental data in Ref. [27]). The experimentally measured J/^ cross section 
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includes the direct J/\l/ component (ctj/^) as well as the decays of higher 
charmonium states Xc and i.e. 



(3) 



Note, we do not distinguish the Xci(l-P) and Xc2(l-P) states. Instead, we use 
only the Xci(l-P) state (which we denote as Xc), however, with an increased 
branching ratio for the decay to J/\l/ in order to include the contribution of 
Xc2(l-P), i-e. B{xc J/^) = 0.54. Furthermore, we adopt B{'^' J/^') = 
0.557 from Ref. [28j. 

We recall that (as in Refs. [T^14f29ll30f31] ) the charm degrees of freedom in 
the HSD approach are treated perturbatively and that initial hard processes 
(such as cc or Drell-Yan production from NN collisions) are 'pre-calculated' 
to achieve a scaling of the inclusive cross section with the number of projectile 
and target nucleons as Ap x At when integrating over impact parameter. For 
fixed impact parameter b, the cc yield then scales with the number of binary 
hard collisions Ncou {cf. Fig. 8 in Ref. t-13j). 

In addition to primary hard A^A^ collisions, the open charm mesons or char- 
monia may also be generated by secondary meson-baryon {niB) reactions. 
Here we include all secondary collisions of mesons with baryons by assuming 
that the open charm cross section (from Section 2 of Ref. [13]) only depends 
on the invariant energy ^/s and not on the explicit meson or baryon state. 
Furthermore, we take into account all interactions of 'formed' mesons - after 
a formation time of rp = 0.8 fm/c (in their rest frame) [32] - with baryons 
or diquarks. For the total charmonium cross sections from meson-baryon (or 
ttN) reactions we use the parametrization (in line with Ref. |33j): 



with 7 = 7.3 and b = 1.24 mb, which describes the existing experimental data 
at low ^/s reasonably well, as seen in Ref. |llj . 

Apart from the total cross sections for charmonia we also need the differential 
distribution of the produced mesons in the transverse momentum pp and the 
rapidity y (or Feynman xp) from each individual collision. We recall that 
Xp = Pz/p^°'^ ~ '^Vzl with "Pz denoting the longitudinal momentum. For 
the differential distribution m. xp from A^A^ and ttA^ collisions we use the 
ansatz from the E672/E706 Collaboration [31] and for the pt distribution a 
power low parametrization from Ref. [35] which has been fixed by the STAR 





(4) 
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Fig. 4. Differential cross section of J /'^ production in pp collisions at ^/s = 200 GeV 
at mid-rapidity {\y\ < 0.35, l.h.s.) and at forward rapidity (1.2 < \y\ < 2.2, r.h.s.) vs 
the transverse momentum squared as implemented in HSD (solid line) compared 
to the PHENIX data from Ref. [23j (dots). 



data [36], i.e. 
dN 



dxpdpj 



(5) 



with = 3.5 GeV/c and Cp^ = —8.3. The exponent c is given by c = 
a/(l + b/y/s) and the parameters a,b are chosen as qnn = 16, b^^ = 24.9 
GeV for A^A^ colhsions and a,rAr = 4.11, b.^^ = 10.2 GeV for vrA^ colhsions. 

The resulting rapidity distribution for J/\l/ production in pp colhsions at ^/s = 
200 GeV is shown in Fig. [3] which is in line with the data from Ref. [23] within 
error bars. We also present the pp —>■ J/'^ + X differential cross section in p"^ 
at mid-rapidity {\y\ < 0.35) and at forward rapidity (averaged in the interval 
1.2 < \y\ < 2.2) in Fig. |H The HSD parametrization is compared to the 
recent measurements of the corresponding quantities by PHENIX [23]. The 
total cross sections for D + D production in this study are the same as those 
presented in Ref. [TT] . 

The parametrizations of the total and differential cross sections for open charm 
mesons from pN and vrA^ collisions are taken as in Refs. [T3lfTi] . apart from a 
readjustment of the parameter a^N in order to reproduce the recently mea- 
sured rapidity distribution of J/\E' in p + p reactions at ^/s = 200 GeV by 
PHENIX [23]. 

For D,D*,D,D* - meson {TT,ri, p,uj) scattering we refer to the calculations 
from Ref. |371i38] which predict elastic cross sections in the range of 10-20 mb 
depending on the size of the form factor employed. As a guideline we use a 
constant cross section of 10 mb for elastic scattering with formed mesons and 
also baryons, although the latter might be even higher for very low relative 
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Fig. 5. J/^' production cross section in d + Au collisions relative to that in p + p 
collisions (see text for the definition of RdA) in HSD (red stars) as compared to the 
PHENIX data [39] (full dots). 

momenta. We will discuss this issue in more detail in Section 15. 1[ 



4 Baryonic ('normal') nuclear absorption 



The yield of J/\l/ in p + A and A + A reactions is modified compared to that in 
p + p scaled with the number of initial binary scatterings Ncou [lOll^ . Indeed, 
the produced cc can be dissociated or absorbed on either the residual nucleus 
of the projectile or target or on light co-moving particles (usually on mesons 
or, at high energy, on partons) produced in the very early phase. The latter 
reactions are only important in nucleus-nucleus collisions and not in p + A 
or d + A a.s the number of 'co movers' created in proton- or deuteron- induced 
processes is small. In contrast, charmonium absorption on baryons is the lead- 
ing suppression mechanism in d + A [p + A) scattering and is an important 
base-line for the study of the absorption in the hot and dense medium created 
in A + A reactions. 

In order to study the effect of charmonium rescattering on projectile/target 
nucleons, we adopt in HSD the following dissociation cross sections of char- 
monia with baryons independent of the energy: 

a,sB = 4.18 mb; (6) 
<yj/^B = 4.18 mb; ct^^b = 4.18 mb; a^,'B = 7.6 mb. 
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In ([6]) the cross section accB stands for a (color dipole) pre-resonance (cc) - 
baryon cross section, since the cc pair produced initially cannot be identified 
with a particular charmonium due to the uncertainty relation in energy and 
time. For the life-time of the pre-resonance cc pair (in it's rest frame) a value 
of Tec = 0.3 fm/c is assumed following Ref. [12]. This time scale corresponds 
to the mass difference of the and J/\l/. 

The values for the cross sections <Jj/i5iN,o'ccN at RHIC energies are currently 
debated in the literature. On one side, all the data on the J/\l/ production in 
p + A at energies ^/s < 40 GeV were found to be consistent with an energy- 
independent cross section of the order of 4 — 7 mb [40j43f44ll45ll46j . On the 
other hand, the corresponding cross sections at the much higher energy of 
a/s = 200 GeV, e.g. at RHIC, are expected to be smaller [17], since part of 
the suppression might be attributed to other (initial-state) cold-matter effects, 
such as gluon shadowing P8ll49|50j . radiative gluon energy loss in the initial 
state or multiple gluon rescattering. We recall that 'shadowing' is a depletion 
of low-momentum partons in a nucleon embedded in a nucleus compared to 
the population in a free nucleon, which leads to a lowering in the charmonium 
production cross section. The reasons for depletion, though, are numerous, 
and models of shadowing vary accordingly. There is, therefore, a considerable 
(about a factor of 3) uncertainty in the amount of shadowing predicted at 
RHIC [i8lli9]l50ll5T]l52] . In the analysis of the d + Au data at = 200 GeV, in 
which the maximum estimate for the effect of the shadowing was made [TTIfSO] . 
the additional absorption on baryons allowed by the data was found to lead to 
(^j/<i!N = 1 — 3 mb or higher, if some contribution of ant i- shadowing is present. 
The authors of [50] advocate (Jj/^n = 3 mb in order to preserve the agreement 
with the data of the Fermilab experiment E866. The PHENIX Collaboration 
[39] finds a breakup cross section of 2.8li;J (using EKS shadowing) which 
still overlaps with the CERN value of 4.18 mb (though with large error bars). 
However, the theoretical uncertainty is still large, since in the works above 
only an approximate model for baryonic absorption was applied and not a 
microscopic transport approach that e.g. also includes secondary production 
channels of charm pairs as described in Section 3. 

Within HSD we have found the baryoninc absorption cross sections ([6]) to 
agree with the data at SPS energies jlT]. In Fig. [5] we compare the HSD result 
(employing the same cross sections (l6|) for baryonic absorption and neglecting 
shadowing) for the J/\l/ production in d+Au collisions at -y/s = 200 GeV to the 
inclusive PHENIX data [39]. The quantity plotted is the nuclear modification 
factor defined as 





Fig. 6. The ratio RdA © for backward, central and forward rapidity bins as a 
function of the number of binary colhsions NcoU foi" d + Au at y/s = 200 GeV. The 
experimental data have been taken from Ref. [39j. The HSD results (stars connected 
by red dashed lines) show calculations without including low-x gluon shadowing and 
slightly overestimate RdA in the forward interval 1.2 < y < 2.2. The theoretical error 
bars are due to the finite statistics of the calculation. 



where dNj^^/dy is the J/"^ invariant yield in d + A collisions, dNY/^/dy 
is the J/\l/ invariant yield in p + p collisions; (Ncoii) is the average number 
of binary collisions for the same rapidity bin. In our analysis we have used 
(N^^ii) = 7.6 ± 0.3 according to the PHENIX estimate [39]. 

It is seen from Fig.[5]that the calculations follow approximately the decrease in 
RdA with rapidity, however, with a tendency to overshoot at forward rapidity. 
Within error bars we find the values of accs from to be compatible with 
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the inclusive RHIC measurement as well as with the lower energy data |44j . 
This finding is also in line with the analysis of the PHENIX Collaboration in 
Ref. [32] 

In order to shed some further light on the role of shadowing, we compare our 
calculations for RdA in different rapidity bins as a function of the centrality of 
the d + Au collision, which in Fig. [6] is represented by the number of binary 
collisions Ncoii- The latter number is directly taken from the number of binary 
hard NN collisions in the transport calculation while the comparison with 
experiment is based on a Glauber model analysis of the data similar to that 
performed in Ref. [S3] • The actual results displayed in Fig. El (stars connected 
by dashed lines) and the PHENIX data from Ref. [39| are roughly compatible 
for the rapidity intervals -2.2 < y < -1.2 and \y\ < 0.35, but demonstrate that 
the suppression at forward rapidity (1.2 < y <2.2) is underestimated in the 
color-dipole dissociation model with a constant cross section of 4.18 mb. This 
clearly points to the presence of shadowing effects at least at forward rapidities 
which is not so pronounced in the inclusive data set in Fig. [51 A more serious 
question is a quantification of the shadowing due to the limited statistics of 
both the experimental data and the calculations. Here we do not attempt to 
attribute a fixed number for the shadowing effect but merely point out that 
independent high statistics data will be necessary to fix this unsatisfactory 
situation from the experimental side. 

Nevertheless, some note of caution is appropriate for the further analysis of 
charmonium suppression in Au + Au collisions: There are 'cold nuclear mat- 
ter effects' such as 'gluon shadowing' beyond those incorporated in the trans- 
port calculations, and especially quantitative statements about any 'agreement 
with data' might have to be reconsidered. In case of Au + Au reactions the 
shadowing from projectile/target will show up symmetrically around y = 
and in part contribute to the stronger J/"^ suppression at forward/backward 
rapidities. Nevertheless, following Granier de Cassagnac [53], an anomalous 
suppression of J/\E' beyond 'cold nuclear matter effects' is clearly present in 
the Au + Au data to be investigated below. 



5 Scenarios for the anomalous J/\I' suppression 

It is well known that the baryonic (normal) absorption alone cannot explain 
the suppression of charmonia in heavy-ion collisions with increasing centrality 
[7j. We have implemented in HSD several different mechanism for the addi- 
tional (anomalous) suppression of charmonia which will be explained in the 
following Subsections. By comparing the results from these scenarios to each 
other and to the available data the mechanism of charmonium interactions 
with the medium can be probed. 
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5.1 ^Comover' suppression (and recombination) 



First of all let us stress that the interactions with 'comoving' mesons lead 
not only to the dissociation of charmonia, but also to their recreation via the 
inverse recombination process D + D — > cc+m, where m = {vr, p, u, K, ...}. As 
already pointed out before, the J/\l',Xc5^' formation cross sections by open 
charm mesons or the inverse 'comover' dissociation cross sections are not well 
known and the significance of these channels is discussed controversely in the 
literature [Hl6j54ll55ll56f57ll58j . We here follow the concept of Refs. p3p5] and 
introduce a simple 2-body transition model with a single parameter |Mop, that 
allows to implement the backward reactions uniquely by employing detailed 
balance for each individual channel. 

Since the charmonium-meson dissociation and backward reactions typically 
occur with low relative momenta ('comovers'), it is legitimate to write the 
cross section for the process 1 + 2— >-3 + 4as 




where Ek denotes the energy of hadron k [k = 1,2,3,4), respectively. The 
initial and final momenta for fixed invariant energy ^/s are given by 



2 _ (s - (mi + m2)^)(s - (mi - m2) 

Pi 



As 

^2 ^ (g - (mg + m4f){s - (mg - m4f) 

where m^ denotes the mass of hadron k. In ([8]) |Mjp (i = Xc J ■, stands 
for the effective matrix element squared, which for the different 2-body chan- 
nels is taken of the form 



\Mi\^ = \Mi\^ for (7r,p) + (cc)i^D + D (10) 
|MiP = 3|Mi|2 for {■K,p) + {cc)i^ D* + D, D + D*, D* + D* 

|M,|2 = i|M,p for {K,K*) + {cc),^D., + D, D., + D 

|M,|' = |M,|2 for {K,K*) + {cc)i^ D, + D*, D^ + D*, Dl + D, 

Dl + D, Dl + D* 

The relative factors of 3 in (fTOj) are guided by the sum rule studies in [59] 
which suggest that the cross section is increased whenever a vector meson D* 
or D* appears in the final channel while another factor of 1/3 is introduced 
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for each s or s quark involved. The factor {{m^ + m^) / y/s) in ([S]) accounts for 
the suppression of binary channels with increasing ^/s and has been fitted to 
the experimental data for the reactions n + N ^ p + N,uj + N,(l) + N, + A 
in Ref. [60]. 

We use the same matrix elements for the dissociation of all charmonium states 
i {i = Xc, J ) with mesons: 



Mj/*r = \M^f = |M^f = \M^\\ (11) 



We note for completeness that in Ref. [H] the parameter |Mop was fixed by 
comparison to the J/\l/ suppression data from the NA38 and NA50 Collabo- 
rations for S+U and Pb+Pb collisions at 200 and 158 AGeV, respectively. In 
a later study |Tl], however, this parameter has been readjusted in accordance 
with the updated value of the cross section ([6]) of charmonium dissociation 
on baryons (following the latest NA50 and NA60 analysis [I0P3] ). The best 
fit is obtained for |Mop = 0.18 fm^/GeV^; this value will be employed in our 
following studies, too. 

The advantage of the model introduced in [H|IT5] is that detailed balance for 
the binary reactions can be employed strictly for each individual channel, i. e. 



(2gi + l)(2g2 + l) Pi 
(2^3 + 1) (2^4+1) pY 



0-3+4^1+2(5) = q-i+2^3+4(-s) ,^c^ , , 1^ ri' (12) 



and the role of the backward reactions ( (cc) j+meson formation hy D + D fla- 
vor exchange) can be explored without introducing any additional parameter 
once I Mop is fixed. In Eq. ( IT^ the quantities Sj denote the spins of the parti- 
cles, while pj and denote the cms momentum squared in the initial and final 
channels, respectively. The uncertainty in the cross sections f|T2l) is of the same 
order of magnitude as that in Lagrangian approaches using e.g. SU{4) flavor 
symmetry [371138] , since the form factors at the vertices are essentially un- 
known [SU]. It should be pointed out that the 'comover' dissociation channels 
for charmonia are described in HSD with the proper individual thresholds for 
each channel in contrast to the more schematic 'comover' absorption model in 
Refs. EMU. 



The regeneration of charmonia by recombination of D {D*) mesons in the 
hadronic phase was first studied by CM. Ko and collaborators in jSSj- The con- 
clusion at that time was that this process is unlikely at RHIC energies [5^58|62j . 
On the other hand, it has been shown within HSD [14j that the contribution 
of the D + D annihilation to the produced J/\l/ at RHIC is considerable. 
Moreover, the equilibrium in the reaction J/\l/ + m ^ DD is reached (i.e. 
the charmonium recreation is comparable with the dissociation by 'comov- 
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ing' mesons). The reason for such differences is that the pioneering study [58] 
within the hadron gas model was confined to J/\l/ reactions with vr's into two 
particular DD channels {D + D* and D* + D*). On the contrary, in Ref. |14j 
the interactions with all mesons into all possible combinations of DD states 
have been taken into account. Note that the p- meson density at RHIC is large 
such that the channel with the most abundant p-meson resonance is domi- 
nant. Furthermore, in Ref. [H] the feed down from Xc and decays has been 
considered. The results of |T3] are in accordance with independent studies in 
Refs. p3|64|65|l5B] . Later work within the HSD approach [12j has supported 
the conclusions of Ref. [H] and stressed the importance for DD annihilation 
in the late (purely hadronic) stages of the collisions. 

5.2 'Threshold melting' 

This scenario is based on the idea of sequential dissociation of charmonia with 
increasing temperature p71l68ll69|l7n] . i.e. of charmonium melting in the QGP 
due to color screening as soon as the fireball temperature reaches the dissocia- 
tion temperatures of (~ 2Tc for J/\l/, ^ Tc for excited states, where Tc stands 
for the critical temperature of the deconfinement phase transition). In the early 
approaches the temperature of the fireball has been estimated using e.g. the 
Bjorken formula ([T]). We modify the standard sequential dissociation model 
in two aspects: (i) the energy density is calculated locally and microscopically 
instead of using schematic estimates (c/. section [2]); (ii) the model incorporates 
a charmonium regeneration mechanism (by DD annihilation processes). 

The 'threshold scenario' for charmonium dissociation now is implemented in 
a straight forward way: whenever the local energy density e{x) is above a 
threshold value Sj, where the index j stands for J/\l/, Xci ^' ■, the charmonium 
is fully dissociated to c + c. The default threshold energy densities adopted 
are 

Ej/^ = 16 GeV/fm3 , e^c = 2 GeV/fm^, and = 2 GeV/fm^. (13) 

The dissociation of charmonia is widely studied using lattice QCD (IQCD) 
[7T|72f73ll74f75] in order to determine the dissociation temperature (or en- 
ergy density) via the maximum entropy method. On the other hand one may 
use potential models - reproducing the charmonium excitation spectrum in 
vacuum - to calculate Mott transition temperatures in a hot medium. Both 
approaches have their limitations and the quantitative agreement between the 
different groups is still unsatisfactory: 

• (A) Potential models employ the static heavy quark-antiquark pair free 
energy - calculated on the lattice - to obtain the charmonium spectral func- 
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tions. This leads to the dissociation temperatures [76] 

• (B) The maximum entropy method is used to relate the Euclidean thermal 
correlators of charmonia - calculated on the lattice - to the corresponding 
spectral functions and yields higher dissociation temperatures [71] 

TmeltiJ/"^) = 1.7-2 T„ TmeltiXc) = 1.1" 1.2 

or [72] 

TmeltiJ/"^) > 1.5 T„ TmeltiXc) = LIT,. 

Our earlier analysis of experimental data at the SPS in the 'threshold melt- 
ing' approach [llj lead us to conclude from the observation of a considerable 
amount of J/\E' in the most central Pb + Pb collisions that the assumption of 
a melting of J/\l/ close to contradicts the data. Therefore, the values f[T3l) 
are applied also in the current study. 

5.3 Discriminating hadronic and partonic phases 

Two more scenarios are implemented in our present HSD simulations that 
are closely related to the 'comover suppression' and the 'threshold melting' 
scenarios outlined in the previous sub-sections. The essential difference is that 
the comoving hadrons (including the Z)-mesons) exist only at energy densi- 
ties below some energy density Ecut, which is a free parameter. We employ 
^cut = ~ 1 GeV/fm^, which is equal to the critical energy density Sc for 
the parton/hadron phase transition. This scenario clearly separates 'formed 
hadrons' from possible pre-hadronic states at higher energy densities. Indeed, 
it is currently not clear whether D- or D*-mesons survive at energy densities 
above Ec but hadronic correlators with the quantum numbers of the hadronic 
states are likely to persist above the phase transition [7^. One may speculate 
that similar correlations (pre-hadrons) survive also in the light quark sector 
above such that 'hadronic comovers' - with modified spectral functions - 
might show up also at energy densities above Ec- 

We recall that the concept of (color neutral) pre-hadrons - explained in more 
detail in Refs. [T6IT8] - has been also used in the hadron electroproduction 
studies off nuclei in Refs. [TBIT7] as well as for high px hadron suppression 
[TS] or jet suppression at RHIC energies [IS]. It has been found that the 
pre-hadron concept works well for hadron attenuation in nuclei at HERMES 
energies [T6p!7] but underestimates the high pt hadron suppression \X8\ as 
well as the jet attenuation at RHIC energies [19]. Nevertheless, the amount of 
attenuation due to such pre-hadronic interactions emerged to be about 50% 
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of the experimentally observed suppression at RHIC such that their effect 
might not simply be discarded. It should be stressed that the concept of pre- 
hadrons refers to the string breaking mechanism as described in Refs. [TBIITB] 
and is independent on the energy density. A detailed study on the space-time 
evolution of pre-hadrons and their formation to hadrons for pp collisions has 
been performed by Gallmeister and Falter in Ref. [78] . 

In line with the investigations in Refs. [T81IT9] we also study J/\E' produc- 
tion and absorption in Au + Au collisions at ^/s = 200 AGeV assuming the 
absorption of charmonia on pre-hadrons as well as their regeneration by pre- 
hadrons. This adds additional interactions of the particles with charm quarks 
(antiquarks) in the very early phase of the nucleus-nucleus collisions as com- 
pared to the default HSD approach. Since these pre-hadronic (color-dipole) 
states represent some new degrees-of-freedom, the interactions of charmed 
states with these objects have to be specified separately. 

For notation we define a pre-hadronic state consisting of a quark- ant iquark 
pair as pre-meson fh and a state consisting of a diquark-quark pair as pre- 
baryon B. The dissociation cross section of a cc color dipole state with a 
pre-baryon is taken to be of the same order as with a formed baryon, 

= 5.8 mb, (14) 



whereas the cross section with a pre-meson follows from the additive quark 
model as [T6IIT7] 

diss '^rr'^^^^ (^ 



Elastic cross sections are taken as 

<c-B = 1-9 mb, = la%. (16) 



Furthermore, elastic interactions of a charm quark (antiquark) are modeled 
by the scattering of an unformed D or meson on pre-hadrons with only 
hght quarks as 

a% = 3.9 mb, = ^a^^. (17) 

In this way we may incorporate in HSD some dynamics of quark- ant iquark 
pairs with a medium that has not yet formed the ordinary hadrons. However, 
it has to be stressed that further explicit partonic degrees of freedom, i.e. 
gluons and their mutual interactions as well as gluon interactions with quarks 
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Fig. 7. The J/^ nuclear modification factor Raa (USD for Au+Au collisions at 
^ys = 200 AGeV as a function of the number of participants Npart in comparison to 
the data from |41j for midrapidity (full circles) and forward rapidity (full triangles). 
The HSD results for the purely hadronic 'comover' scenario are displayed in terms 
of the lower (green solid) line with open circles for midrapidity J/^'s {\y\ < 0.35) 
and in terms of the upper (red dashed) line with open triangles for forward rapidity 
(1.2 < \y\ < 2.2). 

and antiquarks, are not taken into account in the present HSD approach. 
Therefore, we do not expect to reproduce any details of the measured J/^ 
yield. The study of this particular model situation is motivated first of all 
by the possibility to assess the conceptual influence of charm scattering on 
pre-hadrons (in the early reaction phase) on the final rapidity distribution of 
the J/\l/'s (see below). 



6 Comparison to data 



In the transport approach we calculate the J/\I' survival probability Sj/^ and 
the nuclear modification factor Raa as 



•J J/* — TTTTi"? tiAA 



N'jr N^oU-dN^'' Idy 

where Njl^ and N'^J^ denote the final number of J /'^ mesons and the number 
of J/\l/'s produced initially by BB reactions, respectively. Note that Nj[^ 
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Fig. 8. The ratio of the nuclear modification factors Raa at mid-rapidity (|y| < 0.35) 
and at forward rapidity (1.2 < \y\ < 2.2) vs centrality in Au + Au cohisions at 
y/s = 200 GeV. The HSD results in the purely hadronic scenario ('comover ab- 
sorption') are displayed in terms of the blue dashed line (with open circles) and in 
case of the 'threshold melting' scenario in terms of the violet dot-dashed line (with 
open squares). The error bars on the theoretical results indicate the statistical un- 
certainty due to the finite number of Monte-Carlo events in the calculations. The 
lower full green dots represent the data of the PHENIX Collaboration ^41j . Note 
that the data have an additional systematic uncertainty of ±14%. The lower solid 
(red) line with stars gives the result for the 'comover absorption' scenario when 
including additional pre-hadronic interactions with charm (see text). 



includes the decays from the final Xc- In fflSl) . dN^^ /dy denotes the final 
yield of J/\l/ in AA collisions, dN^^'^ /dy is the yield in elementary pp reactions 
while Ncoii is the number of initial binary collisions. 



The suppression of charmonia by the 'comover' dissociation channels within 
the model described in [TT| for a matrix element squared |Mop = 0.18 fm^/ GeV^ 
has been presented already in Ref. [12j as well as the results for the 'threshold 
melting scenario' employing the thresholds ej/ij, = 16 GeV/fm^, = e^,' = 2 
GeV/fm^. Note that the charmonium reformation channels hj D + D channels 
had been incorporated, too (c/. Ref. [H]). Since the PHENIX Collaboration 
has released a new data set we compare our calculations with the most recent 
PHENIX data [41j in Fig. for the J/^f nuclear modification factor Raa (HHD 
for Au-|-Au collisions at ^/s = 200 AGeV as a function of the number of partic- 
ipants Npart for midrapidity (full circles) and forward rapidity (full triangles). 
The HSD results for the purely hadronic 'comover' scenario are displayed in 
terms of the lower (blue solid) line with open circles for midrapidity J/'^'s 
{\y\ < 0.35) and in terms of the upper (red dashed) line with open triangles 
for forward rapidity (1.2 < \y\ < 2.2). The numerical results appear accept- 
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Fig. 9. The J/^' nuclear modification factor Raa (|18P for Au+Au collisions at 
^/s = 200 AGeV as a function of the number of participants Npart in comparison 
to the data from [51] for midrapidity (full circles) and forward rapidity (full trian- 
gles). The HSD results for the hadronic 'comover' scenario including additionally 
pre-hadronic interactions of charm according to (I14p - ()17p are displayed in terms of 
the upper (green solid) line with open circles for midrapidity J/'if's {\y\ < 0.35) and 
in terms of the lower (orange dashed) line with open triangles for forward rapidity 
(1.2 < \y\ < 2.2). 

able at midrapidity {\y\ < 0.35) but the even larger suppression at forward 
rapidity (seen experimentally) is fully missed (cf. Ref. [T2]). 

The failure of the traditional 'comover absorption' model as well as 'threshold 
melting' scenario at the top RHIC energy is most clearly seen in the centrality 
dependence of the ratio of the nuclear modification factors Raa at forward 
rapidity (1.2 < \y\ < 2.2) and at mid-rapidity {\y\ < 0.35) as shown in Fig. [HI 
The HSD results in the purely hadronic scenario ('comover absorption') are 
displayed in terms of the blue dashed line (with open circles) and in case of 
the 'threshold melting' scenario in terms of the dot-dashed violet line (with 
open squares). The error bars on the theoretical results indicate the statistical 
uncertainty due to the finite number of Monte-Carlo events in the calcula- 
tions. The lower full green dots in Fig. [8] represent the corresponding data of 
the PHENIX Collaboration ^Tj which show a fully different pattern as a func- 
tion of centrality (here given in terms of the number of participants Npart)- 
The failure of these 'standard' suppression models at RHIC has lead to the 
conclusion in Ref. [12] that the hadronic 'comover absorption and recombina- 
tion' model is falsified by the PHENIX data and that strong interactions in 
the pre-hadronic (or partonic) phase should be necessary to explain the large 
suppression at forward rapidities. 
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Fig. 10. The rapidity distribution dNj/^/dy for different centralities from the stan- 
dard 'comover' model (dashed blue lines) and the 'comover' model with additional 
pre-hadronic interactions of charm according to (fH|) - (fT7|) (solid red lines) . The full 
dots show the respective data from the PHENIX Collaboration [H]. The calculated 
lines have been smoothed by a spline algorithm. The reactions are Au+Au at ^/s 
= 200 GeV. 

In this work we follow up the latter idea and incorporate in the 'comover 
scenario' the additional pre-hadronic cross sections (fT^ - f lT7|) for the early 
charm interactions to have a first glance at the dominant effects. The J/\& 
suppression pattern in this case is shown in Fig. [H] in comparison to the same 
data as in Fig. [71 Now, indeed, the suppression pattern for central and forward 
rapidities becomes rather similar to the data within the statistical accuracy of 
the calculations. Indeed, the ratio of Raa at forward rapidity to midrapidity 
now follows closely the experimental trend as seen in Fig. [8] by the lower red 
solid line. 

Some further information may be gained from the J/"^ rapidity distributions 
in Au+Au collisions at RHIC. The latter distribution is shown in Fig. [10] in 
comparison to the PHENIX data for central collisions (upper l.h.s.), semi- 
central (upper r.h.s.), semi-peripheral (lower l.h.s.) and peripheral reactions 
(lower r.h.s.) for the standard 'comover' scenario (dashed blue lines) and the 
'comover' model including additionally pre-hadronic interactions of charm ac- 
cording to (fT^ - (fTTjl (solid red lines). Whereas for peripheral reactions these 
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additional early interactions practically play no role, the latter lead to a nar- 
rowing of the J/\Ef rapidity distribution with the centrality of the collision 
(roughly in line with the data). In the standard 'comover' model an opposite 
trend is seen: here the interactions of charmonia with formed hadrons produce 
a dip in the rapidity distribution at ?/ ~ which increases with centrality since 
the density of formed hadrons increases accordingly around midrapidity. Since 
the total number of produced cc pairs is the same (for the respective central- 
ity class) and detailed balance is incorporated in the reaction rates we find an 
surplus of J /'^ at more forward rapidities. The net result is a broadening of 
the J/"^ rapidity distribution with centrality opposite to the trend observed 
in experiment. 

Summarizing the results displayed in Figs. 7 - 10 we like to point out that 
the hadronic 'comover' dynamics for charmonium dissociation and recreation 
- as well as the standard charmonium 'melting' scenario - do not match the 
general dependences of the J/\E' in rapidity and centrality as seen by the 
PHENIX Collaboration. In fact, a narrowing of the J /'^ rapidity distribution 
cannot be achieved by comover interactions with formed hadrons since the 
latter appear too late in the collision dynamics. Only when including early pre- 
hadronic interactions with charm a dynamical narrowing of the charmonium 
rapidity distribution with centrality can be achieved as demonstrated more 
schematically within our pre-hadronic interaction model. Consequently, the 
PHENIX data on J/\I' suppression demonstrate the presence and important 
impact of pre-hadronic or partonic interactions in the early charm dynamics. 
This finding is line with earlier studies in Refs. [H|18|19j demonstrating the 
necessity of non-hadronic degrees of freedom in the early reaction phase for 
the elliptic flow f2, the suppression of hadrons at high transverse momentum 
Pt and far-side jet suppression in central Au+Au collisions at RHIC energies. 



7 Predictions, excitation functions and comparison to the statisti- 
cal hadronization model 

In this Section we continue with predictions for future measurements as well 
as model comparisons in order to allow for an experimental discrimination 
between the model concepts. 

7.1 ^' as an independent probe 

As pointed in Ref. an independent measurement of \E'' will provide further 
information on the charm reaction dynamics and final charmonium formation. 
For instance, a leveling off of the ^' to J/\E' ratio with increasing centrality 
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Fig. 11. Survival probability of ^' in Au + Au reactions at ^/s = 200 GeV in the 
'threshold melting' (dashed line with triangles) and 'comover' suppression (solid line 
with circles) approaches, see text for details. 



would be a signal for charm chemical equilibration in the medium [31IHI5] . Ad- 
ditionally, it provides a very clear distinction between the 'threshold melting' 
scenario and the 'comover' approach. Since detailed predictions for the to 
J/\l/ ratio as a function of centrality have already been presented in Ref. [11] 
for FAIR and SPS energies we here complement the latter studies by results 
for the top RHIC energy although the suppression of mesons has not yet 
been measured at RHIC. 



In Fig [TT], we accordingly present the survival probability S^^i defined as 



(19) 



for Au + Au at y/s = 200 GeV. In equation (fT9l) . Nj-^^ and A^*^ denote the 
number of final \[^' mesons and of those produced initially by BB reactions, 
respectively. One can see from Fig [11] that the 'threshold melting' scenario at 
RHIC predicts an almost complete melting of while a hadronic 'comover' 
absorption scenario shows a gradual decrease of the number of \I'' with Npart- 
Similar differences between the models have also been found at SPS energies 
[11] where the presently available data sets clearly favor the 'comover' model. 
On the other hand our predictions for the top RHIC energy imply that the 
signal will be very low for mid-central and central Au-|-Au collisions such 
that actual measurements will turn out to be very demanding. 
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Fig. 12. Ratio of the averaged J/^ to vr multiplicity for Pb + Pb at the SPS beam 
energy of 158 A-GeV at mid-rapidity (l.h.s.) and in full 47r acceptance (r.h.s.) as a 
function of the number of binary collisions Ncoii for the different suppression scenar- 
ios implemented in HSD - the 'comover' model (dashed blue line with open circles) 
and the 'threshold melting' scenario (green dot-dashed line with open triangles) - in 
comparison to the statistical model by Gorenstein and Gazdzicki p] (r.h.s.; straight 
orange line) and the statistical hadronization model by Andronic et al. [3J (l.h.s.; 
solid black line). 



40 



S 30 



A 
V 
A 



20 



10 



Au+Au, s =200 GeV, mid-rapidity 




'~~r u * 

•- 0— -o— •° 



Q~--ti 4 ^ 

HSD 

— o— Comover 
—A— QGP threshold 
- Andronic etal — * — Prehadron Interactions 



20 



15 



10 



Au+Au, s'^'=200 GeV, An | - 




^ Gorenstein-Gazdzicki 

/ 

X — o — Comover 

X —A— QGP threshold 

Prehadron interactions 



200 



400 



600 

N , 



800 



1000 1200 



200 400 600 800 1000 1200 



Fig. 13. Same as Fig. [I2] but for Au^ Au at the top RHIC energy of ^ = 200 GeV. 
The red solid line shows additionally the result of the 'comover' model including 
the pre-hadronic charm interactions (see text). 



1.2 Testing the assumption of statistical hadronization 



The assumption of statistical hadronization - i.e. of J/\I/'s being dominantly 
produced at hadronization in a purely statistically fashion according to avail- 
able phase space and the number of available c and c quarks - leads to a 
scaling of the (J/\E')/(/i) ratio with the system size [2], where (h) is the aver- 
age hadron multiplicity. Since {h) ~ (vr), we calculate the ratio {J/'^)/{tt) in 
HSD in the different scenarios for charmonium suppression: 

• 'threshold melting' + recombination via DD ^ cc + m including the back- 
ward reactions cc + m —>■ DD, 
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• hadronic ('comover') absorption: DD cc + m and the backward reactions 
cc + m — >■ DD; 

• 'prehadron interactions': DD cc+m and the backward reactions cc+m 
DD as well as early pre-hadronic charm interactions as described in Section 
6. 

The results of our calculations are shown in Fig. [12] together with the predic- 
tion of the statistical model of Gorenstein and Gazdzicki [2j for the full phase 
space (straight orange line; r.h.s.) and the statistical hadronization model by 
Andronic et al. [3|79] for mid-rapidity (solid black line; l.h.s.) for Pb+Pb at 
158 A-GeV. The centrality dependence here is given by the number of initial 
binary collisions Ncoii- The actual comparison in Fig. [12] indicates that the sta- 
tistical model by Andronic et al. [3J predicts a sizeably larger J /'^ to ti ratio at 
midrapidity for peripheral and semi-peripheral reactions than the microscopic 
HSD results for the different scenarios. For central reactions - where an ap- 
proximate equilibrium is achieved - all scenarios give roughly the same ratio. 
In full 47r phase space the HSD results indicate also a slightly higher J/\E' to vr 
ratio in the 'comover' model relative to the 'melting' scenario but both ratios 
only weakly depend on centrality roughly in line with the statistical model 
of Gorenstein and Gazdzicki [2] (orange straight line). Consequently, only 
peripheral reactions of heavy nuclei might be used to disentangle the different 
scenarios at top SPS energies at midrapidity (or in full phase space). 

The situation is different for Au+Au collisions at the top RHIC energy as 
may be extracted from Fig. [13] where the J/\E' to pion ratio (l.h.s.: at midra- 
pidity; r.h.s.: for 47r acceptance) is shown again as a function of Ncoii- The 
standard 'comover' model (dashed blue lines) is only shown for reference but 
is unrealistic according to the analysis in Section 6. We find that the 'comover' 
model with early pre-hadronic charm interactions (solid red line with stars, 
l.h.s.) is very close to the statistical hadronization model [3] (solid black line) 
at midrapidity except for very peripheral collisions. The 'threshold melting' 
scenario follows the trend in centrality but is down by about 30%. Thus at 
midrapidity there is no essential extra potential in differentiating the scenar- 
ios. Considering the full 47r acceptance (r.h.s.) we find a practically constant 
J/^ to pion ratio for Ncou > 200 from the HSD calculations as expected from 
the statistical model, however, the early model of Gorenstein and Gazdzicki 
[2] is down by about a factor of ~ 10 (and may be ruled out by present data). 



7. 3 Excitation functions 

In this Subsection we present the excitation functions for the J/\E' survival 
probability in Au + Au collisions from FAIR to top RHIC energies in the dif- 
ferent scenarios in order to allow for a further distinction between the different 
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Fig. 14. upper part: The excitation function for the J/^ survival probabihty in the 
'QGP threshold melting + hadronic recombination' scenario (dashed green lines 
with triangles) and the 'comover absorption + recombination' model (solid red 
lines with circles) for central (l.h.s.) and minimum bias Au+Au reactions (r.h.s.) as 
a function of the beam energy. Lower part: The ^' to J/^ ratio for the same reac- 
tions as in the upper part of the figure in the 'QGP threshold melting + hadronic 
recombination' scenario (dashed green lines with triangles) and the 'comover ab- 
sorption + recombination' model (upper solid red lines with circles ). 



concepts. The results of our HSD calculations are presented in the upper part 
of Fig. [H] for the 'QGP threshold melting -|- hadronic recombination' scenario 
(dashed green lines with open triangles) and the 'comover absorption + re- 
combination' model (solid red lines with open circles) for central (l.h.s.) and 
minimum bias (r.h.s.) Au+Au reactions as a function of the beam energy. We 
find that from FAIR energies of 20 - 40 A- GeV up to top SPS energies of 158 A- 
GeV there is no significant difference for the J/ survival probability in case of 
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Fig. 15. Elliptic flow V2 of J/^'s produced in central and peripheral In+In collisions 
at 158 A-GeV beam energy in the hadronic 'comover' mode of HSD (open circle 
and open triangle) compared to the NA60 data [80] represented by black diamonds. 



central collisions. The differences here show mainly up in the full RHIC energy 
range where the 'QGP threshold melting + hadronic recombination' scenario 
leads to a substantially lower J/\l/ survival probabilities. In case of minimum 
bias collisions the 'comover absorption + recombination' model (solid lines) 
leads to a roughly energy independent J/\l/ survival probability whereas the 
'QGP threshold melting + hadronic recombination' scenario shows lower J/^ 
survival probabilities (lower dashed green lines) for laboratory energies above 
~ 100 A- GeV due to a larger initial melting of J/\l/ at high energy density. 



A clearer distinction between the different concepts is offered by the excitation 
functions for the \E'' to J/\I' ratio in Au + Au collisions. The calculated re- 
sults are shown in the lower part of Fig. [H] for the 'QGP threshold melting + 
hadronic recombination' scenario (dashed green lines with open triangles) and 
the 'comover absorption + recombination' model (solid red lines) for central 
(l.h.s.) and minimum bias reactions (r.h.s.). Here the is already melting 
away in central Au+Au reactions in the 'QGP threshold melting' scenario 
at bombarding energies above 40 A- GeV whereas a substantial amount of 
survives in the 'comover absorption + recombination' model. Thus measure- 
ments of suppression at the lower SPS or top FAIR energies will clearly 
distinguish between the different model concepts. 
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l.J^ Elliptic flow of charm 



The elliptic flow of particles defined as 

V2{y,PT) = (^^) (20) 

(with = + Py) provides additional information on the collective cur- 
rents and pressure evolution in the early phase of the complex reaction [81] 
since it is driven by different pressure gradients in case of nonvanishing spatial 

2 2 

anisotropy €2 =< p^f^- >• Since €2 decreases fast during the expansion of a 
noncentral reaction the magnitude of V2 gives information about the interac- 
tion strength or interaction rate of the early medium. 

In Fig. [15] we test the HSD result for f2(J/\t') at SPS in the purely hadronic 
'comover' scenario in comparison to the data for f 2 of the NA60 collaboration 
for In+In collisions [80]. In central collisions the elliptic flow is practically 
zero both in the calculation as well as in the experiment whereas in peripheral 
reactions a nonzero flow emerges. The agreement (within error bars) between 
the theory and the data indicates that in line with the reproduction of the J /'^ 
suppression data [TT] the low amount of f 2 does not point towards additional 
strong partonic interactions. Consequently, the present measurements of J/\E' 
elliptic flow at SPS energies do not provide further constraints on the model 
assumptions. 

The situation, however, is different for the collective flow of D-mesons at top 
RHIC energies. In Fig. [12] we show the elliptic flow of D-mesons produced in 
Au + Au collisions at y/s = 200 GeV as a function of the transverse momen- 
tum Pt in HSD (solid blue line with open circles) compared to the PHENIX 
data [82] on f 2 of non-photonic electrons. Here the elliptic flow of D-mesons is 
clearly underestimated in the standard HSD model (cf. Ref. [35]). Only when 
including pre-hadronic charm interactions - as described in Section 6 - the el- 
liptic flow increases (red line with open stars) but still stays clearly below the 
PHENIX data for p^ < 2 GeV/c. We thus have to conclude that the modeling 
of charm interactions by pre-hadronic interactions - as described in Section 6 - 
does not provide enough interaction strength in the early phase of the collision. 
Quite remarkably this finding is again fully in line with the underestimation 
of high Pt hadron suppression [18j as well as far-side jet suppression [I9l in 
the pre-hadronic interaction model. Independently, also the charm collective 
flow points towards strong partonic interactions in the early reaction phase 
beyond the pre-hadronic scattering incorporated so far. 

Since a large fraction of J/\E''s in central Au+Au collisions at RHIC are created 
hj D — D recombination, the elliptic flow of J/^'s obtained from HSD in the 
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Fig. 16. Elliptic flow of L'-mesons produced in Au + Au collisions at ^/s = 200 GeV 
as a function of pT from HSD (solid blue line with open circles) in comparison to 
the PHENIX data [82] on V2 of non-photonic electrons. The red line with open 
stars shows the HSD result for the V2 of D-mesons when including additionally 
pre-hadronic charm interactions as described in Section 6. 

comover (purely hadronic) case is comparatively small, too, and should not be 
in accord with future experimental data. We consequently discard an explicit 
representation of the J/\l/ elliptic flow at RHIC energies since the calculations 
show the f 2 of charmonium to be very close to the D-meson flow within error 
bars. 



8 Summary 

Our present study essentially completes the investigations of charm produc- 
tion, propagation and chemical reactions within the HSD transport approach 
initiated more than a decade ago [30|31j . The present systematic investigation 
extends earlier work to RHIC energies and clearly shows - as advocated before 
[T2] - that the traditional concepts of 'charmonium melting' in a QGP state 
as well as the hadronic 'comover absorption and recreation model' are in se- 
vere conflict with the data from the PHENIX Collaboration at RHIC energies 
whereas both model assumptions work reasonably well at top SPS energies 

The essential new result of this work is that (at top RHIC energies) we find 
evidence for strong interactions of charm with the pre-hadronic medium from 
comparison to recent data from the PHENIX Collaboration ^41j. In partic- 
ular, pre-hadronic interactions (of unformed hadrons) with charm lead to 
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dramatically different rapidity distributions for J/\l/'s and consequently to a 
substantially modified ratio R^^^'^^'^ {J / to R2^^{J/'^) compared to earlier 
calculations / predictions. 



Further results of the present microscopic transport study may be stated as 
follows: 



• The J/\E' suppression in d + Au collisions at = 200 GeV is only roughly 
compatible with the charmonium absorption on nuclei as observed at SPS 
energies in p+A reactions. We find a clear indication for shadowing effects at 
forward rapidity, but a conclusive answer about the size of this effect is not 
possible due to the statistical error bars in both the experimental data and 
the calculations. A proper answer can only be given by future high statistics 
data that allow to fix the scale of shadowing in a model independent way. 

• The to J/\l/ ratio is found to be crucial in disentangling the different 
charmonium absorption scenarios. This result essentially emerges from the 
early dissociation of above the critical energy density ec ~ 1 GeV/fm^ in 
the 'QGP melting scenario' whereas the in the 'comover model' survives 
to higher energy densities. 

• A comparison of the transport calculations to the statistical model of Goren- 
stein and Gazdzicki [2j (in 47r acceptance) or the statistical hadronization 
model of Andronic et al. [3J (at midrapidity) shows differences in the energy 
as well as centrality dependence of the J/\l/ to pion ratio, which might be 
exploited experimentally to discriminate the different concepts. 

• The collective flow of charm in the HSD transport appears compatible with 
the data at SPS energies, but the data are substantially underestimated 
at top RHIC energies (cf. Fig. 16). This not only holds for the standard 
hadronic comover scenario, but also when including interactions of charm 
with pre-hadronic states (unformed hadrons). Consequently the large el- 
liptic flow V2 of charm seen experimentally has to be attributed to early 
interactions of non-hadronic degrees of freedom. 



The open problem - and future challenge - is to incorporate explicit partonic 
degrees of freedom in the description of relativistic nucleus-nucleus collisions 
and their transition to hadronic states in a microscopic transport approach. 
On the experimental side, further differential spectra of charmonia and open 
charm mesons then will constrain the transport properties of charm in the 
early non-hadronic phase of nucleus-nucleus collisions at RHIC (and possibly 
at SPS or even FAIR energies). 
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